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We develop a random-matrix model of two-dimensional dielectric resonators which combines in- 
ternal wave chaos with the deterministic Fresnel laws for reflection and refraction at the interfaces. 
The model is used to investigate the statistics of the laser threshold and line width (lifetime and 
Petermann factor of the resonances) when the resonator is filled with an active medium. The laser 
threshold decreases for increasing refractive index n and is smaller for TM polarization than for 
TE polarization, but is almost independent of the number of out-coupling modes N. The Peter- 
mann factor in the line width of the longest-living resonance also decreases for increasing n and 
scales as \J~N, but is less sensitive to polarization. For resonances of intermediate lifetime, the 
Petermann factor scales linearly with N. These qualitative parametric dependencies are consistent 
with the random-matrix theory of resonators with small openings. However, for a small refractive 
index where the resonators are very open, the details of the statistics become non-universal. This 
is demonstrated by comparison with a particular dynamical model. 



I. INTRODUCTION 

Two-dimensional dielectric microresonators attract 
considerable attention because of their prospective ap- 
plication as microlasers and single- photon cavities, as 
well as sensors in chemical and biological systems (T], 
0) H, 0, IE @- These systems also allow the study of 
generic properties of partially confined waves in a clean 
setting, and hence provide a route to acquire knowledge 
which can be transferred, e.g., to mesoscopic electronic 
devices ■ The analogy is most complete when compar- 
ing non-interacting electrons in cooled, patterned two- 
dimensional semiconductor devices (which usually leak 
through quantum-point contacts) with passive optical (or 
micro- wave) resonators equipped with small openings. 
For such geometries, semiclassical methods and random- 
matrix theory have provided a deep theoretical under- 
standing of the spectral and transport properties, which 
have been found to agree with numerous experiments and 
numerical computations [1, H, EH, EH G3- In par- 
ticular, it is now well established that the predictions of 
random-matrix theory, designed for wave-chaotic systems 
with strong mode-mixing, are of universal applicability 
when the width of the openings W is much less than the 
linear system size L, implying a long mean lifetime in the 
system. For two-dimensional ballistic geometries, this re- 
quirement can be quantified by comparing the numbers 
N oc Wk of outcoupled channels at wavenumber k to 
the number of channels M oc Lk which are mixed by the 
scattering at the boundaries. Standard random-matrix 
universality then requires M 3> N . 

Dielectric microresonators, however, leak everywhere 
around the interface and are far more open than the 
resonators considered in standard random-matrix the- 
ory. The photonic confinement relies on internal re- 
flection, which only becomes perfect for angles of in- 
cidence above a critical value. As a consequence, the 



effective 'openness' increases with the system size, re- 
sulting in N oc M. Universality is no longer guaran- 
teed, and a more detailed modelling is required. This 
modelling also has to account for the possible presence 
of an active medium, essential for microlasers, which 
have been manufactured in many forms and materials 
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In this paper we develop a quantum-dynamical descrip- 
tion of two-dimensional dielectric microresonators which 
combines wave-chaotic propagation of photons inside the 
resonator with the Fresnel laws for reflection and refrac- 
tion at the interface. The general construction is based 
on a variant of the quantum surface-of-section method 
[H, [H [H [H HI 111. The propagation inside the 
resonator is expressed in terms of an internal scatter- 
ing matrix, which we specify by either using random- 
matrix theory [ToL |34| or a quantum-dynamical paradigm 
of chaotic wave propagation, the quantum kicked rota- 
tor [35l | . Similar models have recently been developed 
for mesoscopic electronic and hybrid-superconducting de- 
vices [36, 37, 38, 39,4(3], which require different boundary 
conditions to open up the system and do not allow for 
amplification. 

The random-matrix and quantum kicked rotator vari- 
ants of the model are explored to investigate the thresh- 
old and the quantum-limited line width of wave-chaotic 
dielectric microlasers. The laser threshold is related to 
the decay rate T of the longest living resonances, and 
can be read off from the imaginary part of the frequency 
(which is complex for the non-hermitian operator describ- 
ing an open system). For a good, conventional laser res- 
onator (with almost-perfect mirrors), the line width (full 
width at half maximum of the Lorentzian line shape) is 
given by the Shawlow-Townes formula [4l| 
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where / is the total output intensity. In an open res- 
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onator, the line width is enhanced with respect to this 
prediction by the so-called Petermann factor K, 



Aw = K /S.ljsti 



(2) 



which can be related to the mutual non-orthogonality 
of the resonance modes obtained from the non-hermitian 
operator [H, [H, [H, HH [H, H3] , and has been studied ex- 
tensively in a wide ra nge of q uantum-optical frameworks 
[S HI |H3, M, HI HI IE! HE HS HS HE HI ■ We investi- 
gate how T and K depend on the refractive index n and 
the size of the resonator (quantified by iV and M), and 
also discriminate between the two possible polarizations 
TE and TM (where the electric or magnetic field lies in- 
side the resonator plane, respectively). The results are 
contrasted with predictions of random-matrix theory for 
resonators with small openings [12, H, EH, ES US EJ • 

We start this paper in Section HT1 with a brief summary 
of the basic concepts involved in the description of two- 
dimensional dielectric microresonators. The quantum- 
dynamical models are formulated in Section |IIH In Sec- 
tion [IV] we present the statistics of decay rates T (de- 
termining the laser threshold) and Petermann factors K 
(determining the line width). These results are obtained 
by numerical sampling of the random-matrix ensemble 
and the parameter space of the quantum-kicked rotator. 
A summary of the results and conclusions is presented in 
Section N\ 



II. TWO-DIMENSIONAL DIELECTRIC 
RESONATORS 

We consider planar dielectric microresonators of a 
small width Az whose material properties are character- 
ized by a refractive index n(uj), which may depend on the 
angular frequency u>. Amplification is modelled within an 
effective medium approach, so that the refractive index 
is complex, with — Im n > proportional to the ampli- 
fication rate. The refractive index is taken as homoge- 
neous within the resonator and unity in the surrounding 
medium. The geometry of the resonator is specified by a 
region D in the (x, y) plane. 

We assume that only the lowest-lying transverse mode 
is excited, which provides the best confinement at the 
planar interfaces and hence results in the longest-living 
resonances. For small wavelengths (A = c/lo <C L, 
where c is the velocity of light and L is the linear dimen- 
sion of the resonator in the plane), the radiation outside 
the resonator is mostly confined to the resonator plane. 
Throughout this plane, the electromagnetic field can then 
be described by a scalar wavefunction ip(r), r = {x,y), 
which represents the electric or magnetic field component 
parallel to the z axis (TM or TE polarization, respec- 
tively). 

The wavefunction ip(x,y) obeys the Helmholtz equa- 
tion 



At the interfaces the wavefunction ip and its normal 
derivative d±ip are continuous for TM polarization, while 
for TE polarization these continuity requirements are ful- 
filled by ip and n~ 2 d±ip. 

For interfaces with no curvature, these boundary con- 
ditions result in Snell's law 



sm?7 = nsmx, 



(4) 



which relates the angle of incidence r\ from a wave ap- 
proaching the interface to the angle of refraction \ of the 
wave within the refractive medium. The reflection prob- 
abilities for both polarizations are given by the Fresnel 
laws (12. 63] 
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From within the medium, total reflection occurs for an- 
gles of incidence x > Xc larger than the critical angle, 
which is determined by the condition sinx c = 1/n. 

For the more complicated interface geometry of a mi- 
croresonator, the radiative properties are encoded in the 
scattering matrix S(u), which relates the incoming and 

outgoing amplitudes ip^ ^ = J2j &ij i>j m a suitably 
chosen, flux- normalized basis of scattering states. The 
poles of the scattering matrix are defined by the condi- 
tion S(oj m ) — 00. For lu — Lu m , the Helmholtz equation 
([3]) permits solutions ipm,( r ) with purely outgoing bound- 
ary conditions, = 0. For a real refractive index, 
these poles lie in the lower part of the complex plane 
Imw m = — r m /2 < 0, where T m is the cold-cavity decay 
rate of the mode. Amplification shifts the poles towards 
the real axis. Within this effective-medium approach, the 
laser threshold is reached when the first pole crosses the 
real axis, so that the purely outgoing wave field of the 
laser can be realized at physical, real frequencies lu = uj m 

mm. 

The active medium within the resonator emits radia- 
tion due to the spontaneous emission of photons and their 
subsequent amplification. Assuming total population in- 
version of the active transition within the medium, the 
frequency-resolved output intensity is given by |66l |67| 
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J(w) = — tr [S\uj)S{lo) - 1] 
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[c 2 V 2 +n 2 (r)w 2 ]i/>(r) = 0. 



(3) 



This expression vanishes for a passive medium (with a 
real refractive index), for which the scattering matrix is 
unitary. 

Close to the laser threshold, the radiation is dominated 
by a Lorentzian peak around lo = Rew m . The width 
of the peak can be associated to the spontaneous emis- 
sion processes, which perturb the amplitude and phase 
of the emitted radiation. The Shawlow-Townes formula 
(HJ) arises when the resonance in the intensity ([6]) is cal- 
culated within Breit-Wigner perturbation theory, which 
is based on the mutually orthogonal modes of the closed 
resonator. Equation @ including the Petermann factor 
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accounts for the non-orthogonality of the modes in the 
open resonator. For loL 3> c, the Petermann factor can 
be expressed in terms of the resonance wave function as 

Bill, 
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Above the laser threshold, the resonator becomes un- 
stable, since the gain in the medium outweighs the losses 
through the interfaces. The steady-state intensity is lim- 
ited by pumping and saturation of the medium, which 
requires stud y o f the non-linear regime as done, e.g., in 
Refs. [5ll I56u6^|. The feedback with the medium stabi- 
lizes the amplitude of the laser, but not the phase, whose 
dynamics gives rise to a finite laser line width. The result- 
ing width is half the width obtained from the cold-cavity 
calculation, which also includes the amplitude fluctua- 
tions [48(. Equations (UJ and j2]) include this reduction 
factor of one-half. 



III. QUANTUM-DYNAMICAL MODEL 
A. Construction 

In order to set up the quantum-dynamical model of a 
dielectric microresonator we separate the motion within 
the cavity, which is assumed to be wave-chaotic, from 
the reflection and refraction processes at the resonator 
interfaces, which is assumed to be governed by the Snell 
and Fresnel laws ^ and (O, respectively. 

The construction follows a variant of the quantum 
surface-of-section method [H [29|, [3(1 HH [HJ which we 
adapt to the specifics of wave propagation in a microres- 
onator, and bares similarities to the methods used in an 
efficient numerical scheme for specific cavities [33j . For 
motivation let us first consider the classical ray dynam- 
ics in the system. The successive internal reflections of 
these rays at the interfaces are conveniently represented 
in terms of BirkhofFs canonical coordinates, given by the 
length along the boundary, q, and the sine of the angle 
of incidence, p = sin^. The internal ray dynamics is 
then reduced to a sequence of points (q n ,Pn)i which are 
generated by an area-preserving map M. : (q n ,Pn) ~ > 
(<Zn+i,Pn+i)- At each encounter with the interface, a 
ray is split into a refracted part, which escapes to the 
exterior, and a reflected part which remains inside the 
resonator. The relative amount of reflection is measured 
by the reflection coefficients (JSJ) , which in Birkhoff coor- 
dinates are written 




\/l — p 2 — \/n 2 — p 2 
\J\ — p 2 + \J n~ 2 — p 2 



n 2 p 2 
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(8a) 
(8b) 



for \p\ < l/n, while R = 1 for \p\ > 1/n. 

In wave optics, a similar separation of the internal dy- 
namics and the encounters with the interface can be car- 
ried out for kL 3> 1. In this limit, the internal wave 
function V (int) = V (inUn) + V (int ' out) can be separated 
into the component -0( int < ln ) which propagates away from 
the interfaces, towards the interior of the resonator, and 
the component -0( mt ' out ) which propagates towards the 
interfaces. Because the separation relies on the propaga- 
tion direction it is best carried out in momentum space, 
corresponding to states with a well defined canonical co- 
ordinate p. The geometry of the collision with the inter- 
face is encoded in the internal scattering matrix F(ui), 
which relates the components of the internal wave func- 
tion, 



(int,out) 



F(u)ip 



(int,in) 



(9) 



This is equivalent to Bogomolny's transfer operator [29j], 
originally introduced for closed systems. The matrix 
F(lo) is unitary and symmetric, the latter property 
arising from the time-reversal symmetry in a dielectric 
medium. The dimension M ~ lnt[Cnuj/irc] of F(to) de- 
pends on the perimeter C of the interface, which is pro- 
portional to the system size, C cx L. 

In the classical limit u)L/c oc M — > 00, F(uS) corre- 
sponds to the map M.. When toL/c 3> 1, the momentum 
p is quasi-continuous, 



Pi 



21- M - 1 
2M ' 



I = 1. 



,M. 



(10) 



At the interfaces the internal wave function is cou- 
pled to the external wave field, which can be decomposed 
in the conventional scattering states ?/;( oxt ) = ij;( cxt > ln ) + 
^,(oxt,out)^ rp^g C0U pi m g j s f general form 

^(int,in) = ^(int.out) +r ^(cxt,i„) ; 
^(cxt.out) = _^(oxt,in) +r ^(mt,out) ; (-q^ 

where 1Z is the reflection matrix and T is the transmission 
matrix, constrained by IZ^IZ + T^T = 1. The linear Eqs. 
© and (fTTj) can be solved to arrive at the scattering 
matrix 



S(uj) = -K + TF{u>) 



1 



1 - KF(lu 



■T. 



(12) 



This general form of the scattering matrix has also been 
encountered in recent investigations of electronic trans- 
port and superconductin g h ybrid structures with small 
ballistic openings (S H3, Si IM ■ 

Our construction of the quantum-dynamical model is 
completed by concrete specifications of the matrices 1Z 
and F(ui). For kL 3> 1, Snell's law dictates that the 
matrices 1Z and T in Eq. (fTTj) become diagonal between 
states which conserve np, while the Fresnel laws ||5J) de- 
liver the values of the diagonal elements, 

TZlm = 5lm\jR{Pl), %ra = $imV 1 ~ R (Pl)i ( 13 ) 
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with pi given by Eq. (|T0|) . 

Concerning the internal dynamics encoded in F(to), 
we now follow two routes — we either employ random- 
matrix theory or a generic wave-chaotic quantum map, 
the quantum kicked rotator. 



1. Random-matrix model 

It is commonly accepted that statistical properties of 
wave systems with chaotic ray dynamics are well de- 
scribed by random-matrix theory. In the present case of a 
dielectric system, which preserves time-reversal symme- 
try, this approach amounts to assuming that F at fixed, 
real to can be represented by a unitary symmetric matrix 
randomly drawn from the Circular Orthogonal Ensemble 
(COE) [34]]. This substitution should best capture the 
statistical features of the longest-living resonances, since 
in a realistic resonator ray chaos is only established after 
a few internal reflections. 

In general, the frequency dependence of F(lu) is com- 
plicated. In the subsequent analysis we will only require 
the local dependence around the poles of the scattering 
matrix, and employ as a simple approximation 

F(lj) = exp{inuL/c)U, U from COE(M), fixed, (14) 

where L is the mean propagation distance between two 
reflections. Sabine's law of room acoustics delivers the 
universal expression L = ttA/C, where A is the area of 
the resonator and C is its perimeter. Our approxima- 
tion entails ^tt)" 1 Im tiF^ = MnL/2cn ~ Anu/2-KC. 
This reproduces the correct expression for the mean den- 
sity of states, while a more elaborate modelling of the u> 
dependence would also describe fluctuations around this 
value. 

The random-matrix model of two-dimensional wave- 
chaotic microresonators is obtained by combining Eqs. 
(fT2")l . (fT3")) . and (HU). In order to simplify the notation we 
introduce the quasi-energy 



9 — nujLjc. 

This delivers the scattering matrix 



S 



-K + TU 



exp(-z0) - KU ' 



(15) 



(16) 



where U is an M x M matrix from the COE, while the 
diagonal matrix TZ is given in Eq. (|13[) . 



2. Quantum kicked rotator model 

The quantum kicked rotator [35| is a quantized version 
of the classical standard map on the torus (q,p) <E [0, l) 2 , 
which consists of a combination of torsions p — ► p + q and 
nonlinear 'kicks' q — > q + f(p), where f(p) — f(p + 1) = 
—dv/dp is periodic and can be represented as the deriva- 
tive of a potential v(p). In order to break all symmetries 



apart from time-reversal symmetry we employ a quan- 
tum kicked rotator corresponding to a torsion followed 
by a kick with 



f(p) = K\ sin(27rp) + K 2 sin(47rp + a) 



(17) 



v(p) = — cos(27rp) + — — cos(47rp + a), (18) 

Z7T 47T 

where the kicking strengths K\, K% and the shift a are 
free parameters. The classical dynamics is known to be 
chaotic for K\ % K 2 > 8. 

The propagator, or Floquet matrix, takes the form 



Un 



1 



,S(w-») a -*jI [«(•&)+«($)] 



(19) 



The quantum kicked rotator model of a wave-chaotic mi- 
croresonator follows when this matrix is introduced into 
the scattering matrix (|16p . 



B. Resonances and the Petermann factor 

We now explore the general structure of the scatter- 
ing matrix in the quantum-dynamical model and derive 
general relations for the resonances and the Petermann 
factor. 

The scattering matrix (| 1 2 j) diverges when the resonant 
denominator 1 — 1ZF(lu) has a vanishing eigenvalue. The 
matrix 1ZF(lu) describes one round trip of a wave which 
propagates through the cavity, and then is internally re- 
flected at the interface. This quantization condition is 
analogous to the scattering-matrix quantization condi- 
tion in a closed systems ffuL [7l| , which is recovered for 

n = i. 

For the specific form (| 16(1 of the scattering matrix, the 
quantization condition takes the form 



Kll^m = CX P (-i6» m )* r 



(20) 



which is of the form of an eigenvalue equation with eigen- 
value exp(— i9 m ). The matrix 1ZU can be interpreted as 
a reduced round-trip operator. Since 1ZU is subunitary, 
KU{KUy = TZT^ < 1, all eigenvalues are submodular, 
I exp(— i9 m )\ < 1; therefore the quasienergies 9 m have a 
negative imaginary part, 7 m = — 2Im# m > 0. 

The subunitarity of 1ZU furthermore implies that the 
right eigenvectors ^ m defined by Eq. (|20|) are not mutu- 
ally orthogonal. However, they form a bi-orthogonal set 
with the left eigenvectors, 



($j|# m ) =5 tm , 



(21) 



where the latter are defined by the adjoint eigenvalue 
problem 



$I1W = & m e X p(-i6 m ). 



(22) 



In the original variables, the eigenphases 9 m can be 
interpreted in two different ways: (1) They are poles 
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9 m = nw rn L/c of the scattering matrix for fixed (possi- 
bly complex) refractive index n but variable, generally 
complex uj m , (2) they deliver the threshold condition 
m = n m u)' m L/c of a given resonant state for real u' m 
and a refractive index n m = n' + in' t ' n for which only the 



real part is fixed. For the scattering matrix (fT6|) . both 
problems are intimately related since only the product 
nuj enters the expressions. The cold-cavity poles o> 



are obtained for a real refractive index n = 
threshold amplification is then given by 



(cold) 
m 

The 



n'lm^ old) /Rew 



(cold) 



(23) 



For realistic dielectric microresonators, this relation is 
indeed well established in the semiclassical limit L 3> A = 
Rew/c [7^ | . In such systems, it is moreover reasonable 
to assume that the width of the amplification window of 
the active medium is much less than the center of this 
window, RetUm- Within the amplification window, the 
threshold values n'^ then mainly depend on the imag- 
inary part Imwm°' d ' = —T m /2, where T m is the cold- 
cavity decay rate. The laser threshold is hence deter- 
mined by the longest-living resonances of the cold cavity, 
which are characterized by a small value of r m . 

We now turn to the radiation emitted close to reso- 
nance, brought about by steering the amplification close 
to the threshold value (|23|) . n" — + Sn, while keeping 

lu = Heuim 01 ^ + Suj real. We combine both deviations 
into the quasienergy deviation 9 ss 9 m + 59. The scatter- 
ing matrix (|16p can then be evaluated by only keeping 
the resonant term in the denominator, 



JBr, 



S « TE/tf, 



,-ise 



1 



-$t j 



(24) 



The frequency-resolved output intensity Eq. ^ takes 
the form 



m ^J_ K _ 4sinh 2 (ImM 



where 



2ttc |exp(zRe£0) - exp(Im<S(9)| 2 ' 



A-=(*| n * TO )($t l * TO ). 



(25) 



(26) 



We now can linearize in 69, which produces the 
Lorentzian lineshape 



where the full width at half maximum is Aw 
2<5n Re Wm° ld V 71 '- The total intensity is given by 



(29) 



This recovers relation ([2]) [excluding the factor of 1/2 due 
to the suppression of amplitude fluctuations by the non- 
linear feedback in the lasing regime]. Consequently, K 
defined in Eq. (f2"6")) has to be identified with the Peter- 
mann factor. 

Expression (j2"B"|) assumes that the left and right eigen- 
vectors are normalized according to the bi-orthogonality 
condition (|2"Tj) . When this condition is dropped, the Pe- 
termann factor takes the more general form 



K = 



- m ^ m | 



(30) 



In the present, time-reversal-symmetric situation 
where U = U T , the right and left eigenvectors defined 
by Eqs. (|2T)|) and (f2"2"| can be chosen (via a suitable nor- 
malization) such that <& = {7*\P*. This relation can be 
further exploited by decomposing U = VV T in terms 
of a unitary matrix V , which is fixed up to transforma- 
tions V — > VO with an arbitrary orthogonal matrix O. 
This decomposition allows us to pass to the symmetrized 
eigenvalue problem 

V T 1lV^ m = exp(-t0 TO )tf m , (31a) 
& m V T 1lV = exp(-i9 m )& m , (31b) 

which is solved by the right and left eigenvectors \f = 
V T ^ and $ = V T $ = respectively. In terms of 
the symmetrized eigenvectors, the Petermann factor (|30[) 
takes the form 



K 



\^rn\ 2 ' 



(32) 



This expression is formally analogous to expression ([7]) 
for the Petermann factor in terms of the resonance wave 
function inside the cavity, but involves the eigenvec- 
tors ^ m of the symmetrized reduced round-trip operator 
V T 1ZV. Furthermore, Eq. (|32"|) is formally equivalent to 
the expression used in earlier random-matrix theories for 
cavities with small openings 0, E El Izl Izl I?! 



« ^_4sinh 2 (Im(? m ) 



2ir Re (56» 2 + Im<5<9 2 ' 



(27) 



IV. NUMERICAL RESULTS 



Since we are mainly concerned with long living reso- 
nances in a large cavity, — Ima> m = T rn /2 <C Rew, we 
further linearize the numerator and obtain in the fre- 
quency domain 



1 r 2 

/(a, „ —K , 

2tt (uj-Reut^r + AcuyA 



(28) 



We now employ the quantum-dynamical model with 
scattering matrix (TIT))) in order to investigate the statisti- 
cal properties of wave-chaotic dielectric microlasers. For 
fixed values of the refractive index n and polarization 
(TE or TM), the matrix U is either chosen as a ran- 
dom representative from the COE or as the quantum 
kicked rotator (jTHJ) . The laser threshold follows from 
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FIG. 1: Probability distribution P(y) of the scaled decay rate 
7 m = — 2Im# m , obtained from the eigenvalue equation (|20|l . 
In the upper panels the refractive index is n = 1.5 while in 
the lower panels n = 3.6. The left panels are computed in 
random-matrix theory, while the right panels are computed 
in the quantum kicked rotator model. The matrix dimension 
is M = 200 (solid curves) and M = 100 (dashed curves). The 
labels TM and TE discriminate the different polarizations. 
Each curve is based on 10 5 realizations of U. 



Eq. I]23p. while the Petermann factor follows from Eq. 
(|3"2")l . In order to present the results for the life times 
and laser threshold we use the scaled imaginary part 
7 m = -2Im0 m = T m Z/c. 

In the following we contrast the case of a relatively 
open resonator with refractive index n = 1.5 (close to 
the value of glass) to a relatively closed resonator with 
n = 3.6 (close to the value of Ga(Al)As). In the random- 
matrix model, the computations are based on 10 5 re- 
alizations of U from the COE with matrix dimension 
M = 100 or M = 200 In the kicked-rotator model, 
the same number of matrices with identical dimensions 
are obtained by varying the parameters a, K\ and K%. 
We also present n-dependent averages, which are based 
on 10 4 realizations of U in each of the models. 

Figure [1] shows the probability distribution function 
P(j) of scaled decay rates for both polarizations and the 
two representative values of the refraction index. The 
results obtained for matrices of dimension M = 100 and 
200 are very close, suggesting a mild dependence of P(7) 
on the wavelength when the resonator size increases. For 
the relatively closed resonator, n = 3.6, the distributions 
found in both models are almost identical. For the rel- 
atively open resonator with n = 1.5, however, the limit- 
ing distribution in the random-matrix model is distinctly 
different from the result in the quantum kicked rota- 
tor model. This demonstrates that universality starts to 
break down as the system becomes more open, due to the 
greater influence of short-time dynamics (this breakdown 
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FIG. 2: Probability distribution P(7o) of the extremal (small- 
est) scaled decay rate 70. The symbols and parameters are 
the same as in Fig. [Tj 



has been studied in more detail for the spectrum and res- 
onance wave functions of ballistic systems in [zl, [zll , 
using semiclassical arguments). On the other hand, the 
characteristic qualitative features of the distributions dis- 
play a robust parameter dependence. For both values of 
the refractive index n and in both models, the TE po- 
larization leads to larger decay rates compared to TM. 
For a given polarization, the decay rates decrease with 
increasing n. Both trends are consistent with the gen- 
eral features of the Fresnel laws ([5]) , which provide better 
confinement for TM polarization, and for large refractive 
indices. 

For each realization, the laser threshold is determined 
by the extremal resonance with the smallest rescaled de- 
cay rate, which we denote by 70. Figure [2] shows the dis- 
tribution P(7o) of the extremal decay rate for the same 
parameters as used in Fig. [1] Again the results from 
both models coincide for the relatively closed resonator, 
while for the relatively open resonator the random-matrix 
model has a distinctively more narrow distribution than 
the quantum kicked rotator model. In all cases, TM po- 
larization yields smaller decay rates than TE polariza- 
tion, which is inherited from the behavior of P(7) in Fig. 
Q] The results also show a trend to larger decay rates 
as the dimension M is increased. This trend can be ex- 
plained by the (hardly visible) sharpening of the small-7 
flank in the distribution function P(7) [Fig.[T] when M is 
increased, which suppresses the tail with very small de- 
cay rates. Such a sharpening has also been demonstrated 
in the random-matrix theory of cavities with small ballis- 
tic openings, where the distribution function eventually 
becomes discontinuous [H, [6(| US [U • 

Figure [3] shows the average value (70) of the extremal 
decay rate as a function of the refractive index n. For 
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FIG. 3: Ensemble average (70} of the scaled extremal decay 
rate 70 as a function of refractive index n, for TM and TE 
polarization (squares and circles, respectively). Solid symbols 
are obtained in random-matrix theory, while open symbols are 
obtained in the quantum kicked rotator model. Each data 
point is based on 10 4 realizations of U with matrix dimension 
M = 200. 
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FIG. 4: Probability distribution P(k) of the scaled Petermann 
factor k — K/M . The symbols and parameters are the same 
as in Fig. [T] 



this quantity, the results from the random-matrix and 
kicked-rotator models agree very well. As expected from 
the increasing confinement, (70) decreases with increas- 
ing refractive index n, and is smaller for TM polarization 
than for TE polarization. 

We now turn to the statistics of the Petermann fac- 
tor. Following the results from other variants of random- 
matrix theory @, |H H| |H H [ll H [ll] it has to 
be expected that the Petermann factor increases with 
increasing decay rate of a resonance, and moreover in- 
creases with increasing number of outcoupling channels. 
In the present class of systems this number is simply pro- 
portional to M . 

Anticipating these trends, we show in Fig. 0] the dis- 
tribution P(k) of rescaled Petermann factors n = K/M. 
The distribution represents the Petermann factor of all 
resonances, without discriminating them by width; pa- 
rameters are the same as in Fig. [1] The collapse of the 
curves demonstrates clearly the linear scaling with M. 
Smaller Petermann factors are observed for TM polar- 
ization and/or increasing refractive index, which is con- 
sistent with the reduction of decay rates [see Fig. [T] . In 
contrast to the distribution of decay rates, the results 
from the random-matrix and kicked-rotator model agree 
well for both refractive indices, indicating that the Pe- 
termann factor has a more universal statistics. 

Figure [5] shows the conditional probability distribution 
P(k\j), of Petermann factors with a given value of the 
rescaled decay rate. We again find an approximate lin- 
ear scaling with M. The dependence on polarization and 
refractive index follows the same trends as in Fig. 2] Uni- 
versality is also present, since the results for both models 
are very similar. 

The line width of the microlaser is determined by the 




FIG. 5: Conditional probability distribution P(k|7) of scaled 
Petermann factors k = K/M at a fixed value 7 « 2.0 of the 
scaled decay rate. Results are shown for both polarizations 
and two matrix dimensions. The refractive index is n = 1.5. 



Petermann factor of the longest-living resonance, which 
we denote by Kq and refer to as the extremal Petermann 
factor. (For a given realization, this is not necessarily the 
smallest Petermann factor among all the resonances.) As 
shown in Fig. [6l Ko does not scale linearly with M, but 
scales instead with \[M. This different parametric de- 
pendence can be traced to the fact that the decay rate 
of the extremal resonance is not fixed [Fig. [2]. The same 
parametric dependence, but as a function of TV, is also 
obtained in the random-matrix theory of cavities with 
small ballistic openings, where N <C M [45|, |4(|. Accord- 
ing to Fig. [51 the details of the statistics for the extremal 
Petermann factor are again non-universal for the rela- 
tively open resonator with n = 1.5 — the distribution in 
the random-matrix model is more narrowly peaked than 
in the kicked-rotator model. 

Similarly to what is observed in Fig. [4j the extremal 
Petermann factor in Fig. is reduced when increasing 
refractive index n. However, Kq displays a weaker po- 
larization dependence. These trends are further under- 
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FIG. 6: (color online) Probability distribution P(«o) of the 
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longest- living resonance. The symbols and parameters are 
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lined in Fig. [JJ which shows the ensemble average (ko) 
of the scaled extremal Petermann factor ko = Kq/^/M. 
This average is significantly larger than the modal value, 
which is due to the long tail in the distribution function. 
The average decreases with increasing refractive index, 
similarly to the averaged extremal decay rate in Fig. [3] 
The results from the random-matrix and kicked-rotator 
models converge for large refractive index. 

In summary, Figs. HH3] demonstrate that the laser 
threshold of a wave-chaotic dielectric microresonator de- 
creases for increasing refractive index n and is smaller for 
TM polarization than for TE polarization. The thresh- 
old increases slowly for increasing number of out-coupling 
modes. Figures [6] and [7j show that the Petermann fac- 
tor of a wave-chaotic dielectric microresonator decreases 
for increasing refractive index n and scales as the square- 
root of the number of out-coupling modes, but is less 
sensitive to polarization. Figure Q] shows that this be- 
havior critically depends on the requirement to determine 
the longest-living resonance, which first reaches the laser 
threshold. 

For all quantities, details of the statistics are model- 
dependent for relatively open resonators, but become 
universal for relatively closed resonators with a large 
refractive index (for very closed resonators the results 
eventually converge to the predictions of conventional 
random-matrix theory [l^, 0, HE HE I60l.l6ll|). 



V. CONCLUSIONS 

The quantum-dynamical model of dielectric microres- 
onators presented in this work addresses geometries 
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FIG. 7: Ensemble average (ko) of the scaled extremal Peter- 
mann factor ko = KqI\[~M as a function of refractive index 
n. The symbols and parameters are the same as in Fig. [3] 



which facilitate well-established wave chaos in the semi- 
classical limit A <C L, where L is the typical resonator 
dimension and A is the wave length. We developed two 
variants of the model, one being based on random matrix 
theory while the other is based on the quantum-kicked 
rotator. In general, we found that relatively closed res- 
onators (with a large refractive index) display universal 
statistics while model-dependence sets in for relatively 
open resonators (with a small refractive index). In the 
latter case the non-universal features even affect long- 
living resonances, which classically correspond to rays 
which enjoy many internal reflections at the dielectric 
interfaces. 

The longest-living resonances determine the laser ac- 
tion of the resonator. We concentrated most of our ef- 
forts on the laser threshold and the quantum-limited line 
width, which should be directly accessible in a suitable 
experiment. We found that among the two competing 
polarizations, TM will usually win the mode competition 
for the most stable resonance. The threshold decreases 
with increasing refractive index. The Petermann factor 
in the line width also decreases with increasing refractive 
index, and furthermore scales oc yj L/X. 

Experimentally, Petermann factors of the laser mode 
have been determined for various geometries in Ref. (5?l 
[H, Hil . Recent experimental progress [12] now makes it 
possible to address the life time and Petermann factors of 
individual resonances via direct means. These methods 
are not restricted to the longest living resonance, and do 
not require an active medium (they, hence, also apply 
to passive or absorbing resonators). Such an experiment 
could serve to validate the different scaling cx L/X of 
Petermann factors for resonances of a fixed decay rate. 

Many microresonators of interest have a geometry 
which facilitates ray chaos o nly in parts of the classical 
phase space @, M, M E3, EE HI M M HI M, M, HI, 
l26l \2l\ . While random-matrix theory is not applicable to 
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these systems, the kicked-rotator model can readily ac- 
count for such situations when the kicking strengths are 
suitably reduced. In general, Eq. (fT2| for the scattering 
matrix provides a vehicle to study a wide range of situ- 
ations, including integrable or disordered systems, by a 
suitable round-trip operator 1ZU. Equation (|23| for the 
laser threshold and Eq. (|32p for the Petermann factor do 



not depend on the specific choice of the round-trip opera- 
tor and can serve as a starting point for further analytical 
considerations. 

This work was supported by the EPSRC and by the 
European Commission, Marie Curie Excellence Grant 
MEXT-CT-2005-023778. 
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